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Abstract. Wc compute higher order asymptotic expansions for the partial 
sums of the sequences of central binomial coefficients and Catalan numbers. 
J2k=o Ck ) an d J2k=o^n- We also obtain closed forms for the polynomials 
■ J2t~o { 2 k) xk ana - J2t=o Ck% k over the field of p elements, where q is a power of 

^ \ the prime p. 



1. Introduction 



n + 1 \ n J V n J V n + 1 



Q ■ For n a natural number, the nth Catalan number is defined by 

! n 1 (ln\ (2n\ ( In 

. C„ 

& ■ 

These numbers have a wealth of combinatorial interpretations and applications; 
several dozens of them are collected in |Sta99l Corollary 6.2.3, Exercises 6.19-6.36]. 
As the definition shows, the Catalan numbers are closely related to the central 
binomial coefficients ( 2 ™). Partial sums of these sequences (A000108 and A006134 



in |Slo09j ). and related sums, have recently attracted some interest, especially with 



O 

(S| ■ regard to the values of some of them modulo powers of a prime [PS06[ ISTa[ ISTb] . 

It may be worth noting that, because of the identity ( 2 fc fc ) = Y^=o (J) , the sum 

q ■ Sfc=o ( 2 fc ) equals the square of the Frobenius norm of the lower triangular matrix 

ON . made of the first n + 1 rows of Pascal's triangle. 

In an expanded preprint version of [STaj (larXiv:0709.1665k /5. 23 September 
2007), the authors included some conjectures on Ylk=o (\) an< ^ related sums. One 
^ ■ of them, Conjecture 5.2, concerns the asymptotic behaviour of two such sums; it 

^ ■ claims that 



k=0 v 7 v k=0 v 

The standard notation "~" indicates that the ratio of the two sides tends to one 
as n tends to infinity. In the first part of this note we prove more precise asymp- 
totics than those in Equation (pQ), by means of generating functions and standard 
methods of asymptotic analysis. We also discuss the variations Ylk=o ( 2 ^) cx k and 
^^ =0 Cfca fc of our partial sums, where a is a complex constant. 
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In the second part of the paper we consider the values of those partial sums 
modulo a fixed prime p, that is to say, in a field of characteristic p. Limits and 
asymptotic expansions are generally meaningless in this discrete setting, but there 
is a sense in which an evaluation of an infinite series, such as 



E 

k=0 



f2k\ 1 

, }a k = —= , for a G C with Id < 1/4, 

\k J vT^4^ ' 



has corresponding finite sums over a range < k < q as its prime characteristic 
analogues, where q is a power of p. We illustrate this point of view by giving 
very simple proofs, based on generating functions, of some known congruences 
of this type. They are special cases of more general congruences which were 
obtained in |PS06t ISTa} ISTbj by means of elementary but substantially more 
intricate arguments. By contrast, we only use simple power series manipulations 
and the congruence (1 + x) Q = 1 + x q (mod p). 



2. ASYMPTOTICS 

Equation ([TJ is not hard to prove using the fact that ( 2 ™) ~ 4 n /y / 7m, which 
follows from Stirling's approximation n! ~ n n e~ n y2wn. Indeed, for positive se- 
quences the conditions a n ~ b n and Y^k=i ~^ °°> imply Y^k=i ~ Ylk=i ^k- 
Consequently, we have Ylk=o Ck) ~ Sfc=o v^rA;- The asymptotic behaviour 
of the latter is easily found. For, on the one hand we have J^^ =0 4 fc / \fk > 
(1/y/n) Ylk=o ^ k = (4 n+1 — l)/(3y/n). On the other hand, for any 1 < m < n 
we have 

n A k i n m ~ 1 An+l Am 

yl<±V4 fc + V4 i; < - + — 

t^Vk'V^f^ ft "3v^ 3' 

fc=0 k=m fc=0 

and taking m = \n — log 4 n\ we find that 

^ Vk~ 3^/n-log^n 3n ' 

Taken together, these inequalities imply the former estimate in Equation (TfJ, and 
the latter can be shown similarly. Now we proceed to computing more precise 
asymptotic estimates using the classical method of Darboux. 

Theorem 1. We have 



2k\ 4 n+1 / 1 59 2425 



71 

^ I k j 

k=0 



Proof. Because ( 2 ™) = (— 4) n ( ^ 2 ), the central binomial coefficients have the gen- 
erating function 



n=0 v 7 n=0 v 7 



y/l-4x 
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see [Wil941 Equation (2.5.11)]. It follows that 



which we conveniently write in the equivalent form 



oo 



(3) E^"(E' 

n=0 x fc=0 



having set z = 4x. This function is analytic on the unit disk and has only one 
singularity on its boundary, at z = 1, which is of algebraic type. Because 

4 4 1 4<A 

4^ = 3 " 1 + ((1 - *)/3) = 3 " 3) " (1 ~ ^' 

the Puiseux expansion at z = 1 of the function in Equation ([3]) reads 



1 4 4 



oo 



j=0 



v^T^ 4-z 3 ' 



J-l/2 



According to Darboux's lemma (see [Wil94t Theorem 5.3.1], |Sze75| Theorem 8.4] 
or |Com74l p. 277]), the asymptotic expansion for the coefficient of z n in the 
generating function is formally obtained by adding up the coefficients of z n in each 
term of its Puiseux expansion at z = 1. Truncating and adding the appropriate 
O term we find that 

w 4- 1 (?) - 1 eh»-' ( n " J ; 1/2 ) - 

fc=o ^ 7 i=o ^ 7 

for any nonnegative integer m. Here we have transformed the binomial coefficients 
involved using the standard formula (~ a ) = (— l) n ( n+ ^ _1 ). 

The former estimate in Equation ([I]) follows by setting m = in Equation (j2J), 
using 4 n ( n_1 ' 2 ) = ( 2 ™) = (2n)\/(n\) 2 and Stirling's formula. Of course, higher 
values of m provide sharper results, as we now exemplify. Note that 4 n ( ra ~ n 1/2 ) = 
( 2 ™), whence 



4" 



n-j- l/2\ (2n\ (-l/2)(-3/2) • • • ((-2j + l)/2) 



n J \n J (n-l/2)(n-3/2)---(n-j + l/2)' 
'2ra\ 

J (2n - l)(2ra/3 - 1) • • • (2n/(2j - 1) - l) 
Plugging this into Equation (J4]) multiplied by 4 n we obtain 



(5) ^ V * / 3 V n ; ^ (2n - l)(2n/3 - 1) • • • (2n/(2j - 1) - l) 
+ 0(4 n n" m - 3/2 ), 
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where the summand is interpreted to take the value 1 when j = 0. For example, 
when m = 3 one may use 

2n\ _ _4™_ f J_ _L_ _J3_ 
n ) ~ ^fwn V 8n + I28n 2 + 1024 n 3 + ^ 



(see [GKP94t Exercise 9.60], for example) and obtain the conclusion as given in 
the statement of the theorem after lengthy calculations. It is certainly best to 
let a symbolic manipulation system such as Maple do the calculations, starting 
directly from the right-hand side of Equation (JSJ) . □ 



Theorem 2. We have 

3ru/7m r 8n ' 384n 2 ' 9216n 



n 4" +1 / 5 475 1225 



k=0 

Proof. The generating function for the Catalan numbers, 



(6) £0,*»=-~ ' X 



2x 

n=0 



can be obtained by integrating Equation ([2]) and adjusting the constant term, 
see |Wil94t Equation (2.5.10)]. Proceeding in a similar fashion as in the proof of 
Theorem [TJ we have 



n=0 x fc=0 

which we rewrite in the equivalent form 



oo / n n. 



Vl^ 



n=0 x fc=0 

where z = 4x. Now 



z(4 - z) 



V ' 3=0 j=0 

and hence the Puiseux expansion at z = 1 of our generating function is 



oo 



« \ ( ^17 = Iep- (-y-'wi - *y - (i - -~>' +i/2 ). 

1 J 3=0 

Darboux's lemma tells us that 

4- EC* = -|E(3- (-3)-) (" ~ J ~ 3/2 ) + 0(n— 5 / 2 ), 

fc=0 j=o ^ ' 

for any nonnegative integer m. Note that the coefficient (^) of z n in (1 — z)- 7 in 
the Puiseux expansion gives no contribution to this estimate, because it vanishes 
as soon as n > m; put differently, those terms of the Puiseux expansion add up 
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to a part of the generating function which is analytic at 1. Similar calculations as 
in the previous case lead to 

^ c _2/2n\y. 3. (_i)j + 3-j 

+ 0(4 n n~ m " 5/2 ). 

A Maple calculation with m = A returns the desired estimate. □ 

As concluding remarks, note that this classical method applies similarly to pro- 
duce asymptotics for the modified sums Y^k=o ak ^\k) anc ^ Y^k=o ak ^k, where a is a 
complex constant. We refrain from computing higher order asymptotics for these 
more general sums and limit ourselves to finding the first order asymptotic expan- 
sion for the former sum (the other being similar). A case distinction is necessary 
according to the value of a. 

When | a | > 1/4, the generating function in 



(7) 

n=0 v fc=0 



1 Aa 



y/1 - z 4a 



which generalizes Equation ([3]), still has z = 1 as its dominant singularity (that is, 
the one of smallest modulus). The case m = of the corresponding modification 
of Equation fl4]) then yields 

k (2k\ (4a) n+1 (n - l/2\ ,„ , /2N > (4a) n+1 u ^, , /2N * 



When |a| < 1/4, the pole at z = Aa becomes the dominant singularity of the 
generating function in Equation (J7J). After bringing this singularity to the unit 
circle by the substitution z = Aax, an application of Darboux's method with 
m = yields 



Vo*ffl = — (l + Ojn- 1 )). 



where the square root in the formula refers to the branch of yl — z (on the domain 
\z\ < 1) which evaluates to 1 when z — 0. 

When | a | = 1/4 and a ^ 1/4, the generating function in Equation (J7J) has 
two singularities of modulus one, and the method of Darboux requires adding 
up the contributions from both singularities. If we are content with a first order 
asymptotic expansion, we have 



Wffl = 1 (l + Qjn- 1 / 2 )). 



where the leading term comes from the pole at z — 1, and the other singularity 
contributes to the error term. 
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Finally, when a = 1/4 the function in Equation fl7|) is simply (1 — z) 3//2 . In 
this case our sum admits a closed form 



j%_n V / i,— n \ ' V / \ 



fc=0 x 7 fc=0 

(see sequence A002457 in |Slo09j ). and hence 

n-l 



\ 4 fe =2vn/7r H H - + 0(n~ 

^ Vfc/ V/ \ 8n 128n 2 1024n 3 v 



k=0 

3. Partial sums modulo a prime 
Let q be a power of a prime p. We will show that 

( 8 ) E(?) s fl) (modp) ' and E^^^^ ( mod ^ 

fe=0 ^ ' ^ ' fc=0 



where (|) is a Legendre symbol, and hence is uniquely determined in this case by 
(|) G {0, ±1} and (|) = a (mod p). We will do that by finding closed expressions 

for the polynomials Yll^o Ck) xk anc ^ ^2l=o CkX k over the field of p elements, after 
which the substitution x — 1 will give the desired conclusion. It will be convenient 
to assume that p is odd. However, when p = 2 congruences (JSJ) follow from the 



power series congruences ^^L (\) xk = 1 (mod 2) and ^2^ =Q C] e x k = J^o^ 2 ' -1 
(mod 2), which are easy to prove directly, starting from the definitions of binomial 
coefficients and Catalan numbers. 

Theorem 3. If q is a power of an odd prime p we have 



3-1 

E 

fc=0 



Q fc V = (l_ 4x )(^l)/2 ( modp ). 



Proof. Recall from Equation ([2]) that Ylh=o Cfc) 3 ^ = — 4x)~ 1//2 , and work in 
the formal power series ring Z[[a;]]. Basic facts about binomial coefficients and 
Fermat's Little Theorem imply that (1 — Ax) q = 1 — (4a;) 9 = 1 — 4x q (modp). 
Therefore, noting that all binomial power series involved have integral coefficients, 
we have 

(1 - 4X)- 1 / 2 = (1 - 4x)("- 1 >/ 2 ((1 - AxfY 112 

= (1 - 4x) (9 " 1)/2 (l - Ax q )~ 1/2 (mod p) 

= (1 - 4a;) (9 - 1)/2 (mod (x q ,p)). 

This last congruence, with respect to the modulus (x q ,p), means that the polyno- 
mial, obtained from the power series at the left-hand side by discarding all terms 
of degree q or higher, is congruent modulo p to the polynomial at the right-hand 
side. The desired conclusion now follows. □ 
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Theorem 4. If q is a power of an odd prime p we have 

k 1 - (1 - Ax)^ 2 , , . 

2 , C k x k = — x q 1 (mod p). 

k=0 

Proof. In this case we need to carry the calculation of the previous proof one step 
further, and obtain 

(1 - 4x)~ 1/2 = (1 - Ax) {q ~ l)/2 + 2a; 9 (mod (x q+1 ,p)), 

whence 

(1 - Ax) 1 ' 2 = (1 - 4x)" 1/2 (l - Ax) = (1 - 4a;) (c ' +1)/2 + 2a; 9 (mod (x 9+1 ,p))- 
Therefore, we have 

2_^C n x = — = — x q (mod (x q ,p)), 

71=0 

which implies the desired conclusion. □ 

Congruences (jHJ), for odd p, follow by evaluating on x = 1 the polynomials of 
Theorems [3] and HI using the fact that (— 3)^ -1 ^ 2 = (|) (mod p). This is easy to 
show either by using Jacobi symbols and Gauss' quadratic reciprocity law, or by 
viewing —3 as the discriminant of the polynomial (x 3 — l)/(x — 1) = x 2 + x + I, 
and then evaluating on —3 the quadratic character of the finite field ¥ q . 

Remark 5. It appears that congruences (jSJ) were first established in |PS06l The- 
orem 1.2 and Corollary 1.3]. Those results have the restriction q = p, but have 
(k+d) anc ^ Ck+d in place of Q fc ) and C^, where d is an integer with < d < q. It 
may be possible to extend our method in order to cover this variation. 

The results of [PS06] include congruences for similar sums where ( 2 fc fc ) and Ck 
are multiplied by fixed powers of k, or divided by k. The results of the former 
type can easily be recovered from Theorems [3] and H] by repeated application of 
the operator xD, that is, differentiation followed by multiplication by x. As an 
example, from Theorem [3] we have 

k I ) x k = xD(l - Axf q - l)/2 = 2x(l - Ax) iq ~ 3)/2 (mod p) 

k=l ^ ' 

for p odd. Assuming p > 3 for simplicity and substituting x = 1 we obtain 

^k(^)^2(-3Y q -^ 2 ^^( q -) (modp) 
fc=0 ^ ' ^ ' 

When q = p this is the case d = of the second assertion of [PS06t Theorem 1.2]. 

Remark 6. According to |STal Corollary 1.1], congruences (|S]) actually hold mod- 
ulo p 2 , with a change for p = 3 in case of the former sum. We do not know 
analogues modulo p 2 of our Theorems [3] and HI However, our approach based on 
those results has the advantage of readily allowing an evaluation of the modified 
sums J2l=o { 2 k) ah an d J2k=o^kC^ k modulo p, where a is any algebraic integer. In 
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the special case where a is an ordinary integer prime to p, the former evaluation 
can also be obtained by solving a linear recurrence given in |STbt Theorem 1.1]. 
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